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The electron spectrum structure in the half-filled Hubbard model is considered in terms of the one- 
■ particle Green's functions within many-electron representation. A simple analytical generalization of 

the single-site Hubbard-III approximation is obtained, which takes into account the Fermi excitations 
£SJ . (Kondo terms). The problem of the metal-insulator transition is investigated. The occurrence of a 

three-peak density-of-states structure including the "Kondo" peak at the Fermi level is discussed. 
A comparison with large-d calculations is performed. 
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I. INTRODUCTION 



The problem of strong correlations in many-electron systems is one of the most important in the solid state theory. 
The simplest model to describe correlation effects is the Hubbard model jlj which includes the on-site Coulomb 
interaction. One of most interesting phenomena is the correlation-driven metal-insulator transition (MIT), which 
takes place in a number of transition metal compounds. A simple description of MIT was given by Hubbard Q who 
started from the atomic-level picture and proposed a simple interpolation self-consistent scheme. 

Since the Hubbard works of 60's, a great progress has been achieved in understanding electronic structure of 
highly-correlated systems. Last time, the role of the Kondo effect has been discussed within the large-ei approach (d 
is space dimensionality) which reduces the original periodic Hubbard model to an effective Anderson impurity model 
lS^j' lUtl- Such an approach (dynamical mean-field theory, DMFT) turned out to be rather successful. The corresponding 
density of states (DOS) has three-peak rather than two-peak structure: an additional "Kondo" quasiparticle resonance 
C*~) . at the Fermi level occurs owing to scattering by the local moment. The spectrum structure in large-<i approaches 
is confirmed by the quantum Monte-Carlo (QMC) calculations (see, e.g., ||) and some spectroscopic experimental 
results. Unfortunately, there exist some difficulties in numerical calculations within QMC and large-d approaches, so 
that one needs often to introduce rather high temperatures to resolve these problems. The three-peak structure is 
, not reproduced by most preceding analytical approaches, in particular, by the single-site Hubbard-III approximation 
pL the reason being in that they do not take into account contributions of Fermi-like excitations in a proper way. 
t-H | Thus these approaches do not describe the Brinkman-Rice effective-mass enhancement which is important from the 
experimental point of view. Recently, an attempt has been made to improve the Hubbard-III approximation by 
calculating corrections owing to correlation effects however, the results remained qualitatively unchanged. 

A detailed analysis of Hubbard-III-like approximations was performed in Refs. [^|J^] within the large- z expansion, z 
being the nearest-neighbor number. In the zero order this approach reduces to the simplest Hubbard-I approximation 
0. General expressions for 1/z-corrections in the limit U — > oo were obtained in Ref. ||. The problem of MIT 
within this approach was treated in Ref. j^]. Unfortunately, only a classical approximation (the large-5 1 limit of the 
s-d model which generalizes the Hubbard model) was considered, and the terms with the one-particle occupation 
O ■ numbers, which just describe the Kondo effect in narrow bands jlCjfl , were neglected. 

In the present paper we present a treatment that is based on the method of equations of motion for the many- 
electron Hubbard operators p| , p!2| and is a much more simple than the large-d approach. In Sec. ||, the decoupling 
scheme with account of the Fermi excitations is developed. In Sec. Ill, we present the results of numerical calculations 
and carry out a comparison with previous works. 



II. THE DECOUPLING SCHEME 

We start from the Hubbard model with the electron concentration n = 1 (the half-filled case). The corresponding 
Hamiltonian reads 

U = ^ik4 CT c k(T + U^n^nn, (1) 

ktr i 
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where = cj^c^, c\ a and cj. are the one-electron operators in the Wannier and quasimomentum representation. 
We pass to the Hubbard X-operators 

Xf = \ia){iPl Xf XT' = 8^X9", ^^ = 1, (2) 

so that 

4 = Y.^ c imxf = X? + aXf-°. (3) 

a, 13 

Then the interaction Hamiltonian takes a diagonal form and we obtain 

n = j2 *k(*-°k + ^-rx^ + ^r 2 ) + u e x T- (4) 

kcr i 

Using (^) we have for the one-electron anticommutator retarded Green's function 

G^{E) = (M<L))e = {{X^\clj E + <t{{X^ 2 \cIJ e . (5) 

The energy E is supposed to be referred to the chemical potential which equals U/2 in our case. We write down the 
system of equation of motion 

E((A\B)) E = ({A,B}) + (([A,H]\B)) E 
for the pair of the Green's functions in the right-hand side of (|5|). Using (|^) we obtain in the non-magnetic case 

(E + U/2){{X^\c{ a )) E = 1(1 + t k {(c ka \clj E ) 

+ + 6(X^)) C(lrT \clj E 

q 

+ E + ^ q -<r X kl q |cL))s, (6) 

q 

a(E U/2){{X^\c{ a )) E = 1(1 + t k (( CkCT |cL» B ) 

+ E*q«( 5 (^kf q CT ) + 5 (^k- q ))c qCT |cL»^ 
q 

- E *q«*k-qCq-* + ™U*k%|cL»E: ( 7 ) 

q 

where 5 A = A — (A) is the fluctuation of the operator. Solving the system (^|), (0) we derive 

GUE)=G° k (E)(l-^ a (E)), (8) 



TUE) = E *q«*(*E-q + ^k- q )Cq CT + ^C q _ CT + (TCq.^^Jc^)}^. (9) 



Here 



g °^ = MeW^ W) = E -^ (10) 

is the Green's function of the Hubbard-I approximation (which plays the role of a mean-filed approximation for our 
problem), Fq(E) being the corresponding inverse locator. The Hubbard-I spectrum contains two correlation subbands 
defined by the poles of (n(t) 



£ki,2 = i(*k ± £k), £k = ^U 2 + tl. 
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The Green's function T^^E) describes fluctuation corrections. The corresponding collective excitations are described 
by spin and charge operators. 



of \^ff-ff Q z — ( y++ _ y — 

q — q ' q — 2 q q 

n+ — Y 20 n z — ^ (Y 22 Y 00 °) 

Further we write down the system of equations for the fluctuation Green's functions and perform the decouplings 
which correspond to the first order in the formal parameter l/z (strictly speaking, this expansion is justified in the case 
of long-range electron hopping). For the half-filled band, we have to take into account particle and hole excitations 
in an equal way. However, decouplings can violate the particle-hole symmetry. To preserve this symmetry, we make 
an identical transformation by taking in (^) the Green's functions with symmetrized operator products, e.g., 

<<* k - q %-l C kff)>£ "> ^(<{*k-™ Cq-ff}l4ff)>£ = l((X^C^+C^X k ^\cij E . 

We also use in the equations of motion the Hamiltonian in the symmetrized form, 

ikCkff C kff -> 7^k(Ckff C kff - C k ffC k(7 ). 

Then we obtain for the transverse spin fluctuation contribution 

(E 2 - U 2 /A - Et^-iX^ c^ a }\clj E = (E + C//2)(i q - f k )(/ q - \)G^{E) - U(t k {X^X^) 

'2 



+(i q - t k )(hxi°°, X$-°] aX^X^GUE) 



-(E + C//2)(/ q - i) + U{\[XZ?, X°~°] 

— q q ~ A k-q A -k+q/- 

Note that, as well as the standard Kondo terms, the "many-electron" terms come from the spin-flip processes, but 
not from longitudinal spin fluctuations. For the "transverse" charge contribution we have 

a(E 2 - U 2 /A - M q )<(i{c q _ CT X^ q }|c kCT )) £ = (E + !7/2)(* q + < k )(/ q - \)G^{E) 



uiMx^x^) + (t q + t k )(aXf CT q x q 



+ 2^-q^q 2 ]))Gkff(^) 

+(^ + [//2)(/ q -i)-C/<i[X^ q ,Xf] 

_CrA -q A q ~ A -k-q A k+q/- 

A symmetry of spin and charge degrees of freedom occurs for a symmetric conduction band. 

Further we neglect q-dependence of spin and charge correlations functions and replace them by single-site averages, 
so that 

{S^S-")=2{S^} = {X^), 

Such an approximation is made (although as a rule implicitly) in practically all works on the MIT problem. This 
corresponds to neglecting dynamics of low-energy Bose excitations and may be justified not only in high-temperature 
limit, but also within the 1/z-expansion. For the local-spin subsystem, this approximation is in spirit of the mean-field 
theory. Main part of charge dynamics (the Hubbard splitting U) is also already taken into account in the zero-order 
(Hubbard-I) approximation. A consistent consideration of dynamics can be made in higher orders in l/z. This may 
lead to a change in details of the MIT picture. However, this problem is rather difficult. 
Taking into account (O) we can use the sum rule in (0) to obtain 
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Gk,E) - KE)-a(E)t t - F(E)-t k - F{E) -W)' ' 

•^-^Iw^^ti^^ws^- (13a) 

We have substituted here the one-particle correlation functions in the Hubbard-I approximation, 

1 

2el 



-q 

-a 



(c^X-^) = _[£ ql /(£ ql ) - £ q2 /(£ q2 )], 

/q = (4,^) = -[(^ql - £V2)/(#ql) - (#q2 " l//2)/(£7 q2 )]. 

Due to the symmetry of the bare band, we have 

£/(£q, 2 )<i>(iq) = - /(^qlPHq)- 

q q 

As a result of our way of decoupling, we have in the sums f(Eqi) f(Eqi) — 1/2, and the DOS of the interacting 
system remains symmetric. 

To obtain the self-consistent (SC) approximation we have to replace in ([!(]) the Hubbard-I Green's functions by 
the exact ones, 

and the Fermi functions / q by the exact occupation numbers n q , according to the spectral representation, 



l q = -~ / dEf(E)ImGq(E). 



Then we have the SC equation for the one-electron Green's function in the form (|12]) with 



O Jj2 e\TT 

W = 1 + Jm E ^ G ^ E ) + T E *q G q(^)"q. ( 14& ) 



q 



or/ 

6(25) - F (E) + — ^ * q G q (£> q . (14b) 
q 

Unlike the simplest self-consistency scheme considered in Ref. S [see eq. (32) of that paper], the approximation 
(|l4|), as well as the standard Hubbard-III approximation, does not result in a violation of analytical properties. 

III. RESULTS AND DISCUSSION 

To investigate the MIT problem, we calculate the single-particle density of states 

N(E) = --l m y G k (£). 

k 

The results for the approximations ([l3]) and ( |l4| ) are shown in Figs. [l]-[|, and the critical values for MIT are given in 
Table [j]. The numerical calculation were performed for the square and two cubic lattices with a symmetric bare DOS. 
We also treat the Bethe lattice, i.e. the model semielliptic bare conduction band with 
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4 / (2E X 2 



{W is the bare bandwidth), the rectangular DOS, and the Gaussian DOS 

N(E) = J_ exp 

which corresponds to the hypercubic lattice in the large-d limit. The Gaussian DOS does not have band edges, so 
that the parameter W is determined from the second DOS moment, 




W = Vpa, M2 = J E 2 N{E)dE. 



Then the quantity [12 has equal expressions in terms of W for the Gaussian and semielliptic bare DOS's. 

It is important that the quantity F(E), unlike Fq(E), does not diverge at E — ► 0, i.e. in the centre of the band. 
This fact is just due to many-electron corrections. Therefore the non-self-consistent (NSC) formulas ([l3]) yield a 
metal-insulator transition at U =/= 0, unlike NSC local approximations. However, the corresponding critical value 
U^ sc is rather small. The "false" singularities at the edges of the Hubbard-I bands occur at large U in our NSC 
approximation (see the discussion in Ref. [||). 

The critical value for MIT U c in the standard Hubbard-III approximation Q for an arbitrary bare DOS is given by 

§ 

Uf = 2^i/T 2 . (15) 

The critical value in the SC approximation ( Jl4] ) is changed somewhat in comparison with the Hubbard-III result (see 
Table |). Unlike Ref. 0, where the critical value was decreased by fluctuations, U c /W = 0.67 for the Bethe lattice, our 
approach yields an opposite tendency, in agreement with the results of the QMC calculations at finite temperatures, 
Uc/W ~ 1 (see Table ||). Within the "linearized" DMFT Q, an analytical expression U c can be obtained, which 
exceeds the Hubbard-III value, 

C/ C L = V3U? = 6^. (16) 

As follows from comparison with Table ||, this approximation seems to overestimate somewhat U c . 

The account of the Fermi excitations results in a modification of the DOS form (cf. Ref. ||). In comparison with 
the Hubbard-III approximation, a pronounced pseudogap exists near MIT at U < U c . The same feature can be seen 
from the results of Ref. 0. At small U, a three-peak structure can be seen in Figs. |l|-^, which becomes smeared 
with approaching MIT (the central peak becomes wide, and a pseudogap occurs). The three-peak structure is more 
pronounced in the NSC approximation. The details of our MIT scenario differ from the DMFT picture where the 
central quasiparticle peak is expected to shrink gradually at U — > U c — 0. Probably, this discrepance is connected with 
the overestimation of the role of the damping in our approach. The correct treatment of the damping is a difficult 
problem. To avoid this problem, most large-d calculations are performed at finite temperatures. It should be noted 
that various DMFT versions give somewhat different MIT pictures. In some calculations, a pseudogap develops in the 
metallic phase near MIT (see, e.g., Refs. p4] , p^[ ), as well as in our picture. The occurrence in some works of a very 
high and narrow peak with lowering temperature, even not too close to U c , is probaly an unphysical drawback. In a 
number of large-d calculations (see Refs. (|,|l^,^,^0|) the height of the peak does not change when approaching MIT. 

In the case of square lattice, the situation is more complicated owing to the Van Hove singularity at the band 
centre. The underestimation of U c in the Hubbard-III approximation is confirmed by the calculations of Ref. pj] 
where U c ~ 1.5IF. More weak Van Hove singularities are present for cubic lattices. One can see from Table 1 that the 
difference between our calculations and Hubbard-III results becomes rather strong for these lattices . 

To conclude, we have demonstrated that a simple decoupling scheme enables one to reproduce the non-trivial 
spectrum structure in the half-filled Hubbard model. Our approach yields a qualitative agreement with the results 
of large-d approaches and QMC calculations. At the same time, this can be easily applied for arbitrary two- and 
three-dimensional lattices. In principle, the many-electron Hubbard operator method enables one to consider in a 
regular way the problem of electron structure of systems with the Hubbard splitting. Various types of slave boson 
and fermion representations combined with diagram techniques can be used to this end. 

Since our approach starts from Hubbard's subbands and includes large incoherent contributions, this does not 
reproduce properly the Fermi-liquid (FL) description of quasiparticle states. An account of low-energy spin and 
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charge dynamics would be useful to describe the electron spectrum picture in more detail. A possibility of the 
transition from FL to non-FL behavior which can take place near MIT should be also taken into account. 
The work is supported in part by the Grant of RFFI No. 00-15-96544 (Support of Scientific Schools). 

TABLE I. Critical values of metal-insulator transition for different bare DOS forms in the Hubbard-III approximation, U^, 



"linearized" DMFT, C/ C L , and NSC and SC approximations ^ and (0), [7 C NSC and C/ C sc . 

DOS Uf/W Ut/W U? sc /W U? c /W 

rectangular 1 1.73 0.99 1.22 

semielliptic \/3/2 = 0.866 1.5 0.87 1.06 

Gaussian x/3/2 = 0.866 1.5 0.87 1.06 

square 0.866 1.5 0.87 1.06 

simple cubic 0.707 1.22 0.76 0.99 

bcc 0.612 1.06 0.67 0.92 



TABLE II. Critical values for the metal-insulator transition for the Bethe lattice, [/<? , and in the large-d case, [/«?, from 
different works. 



U C B /W 



Refs. 



Method 



1.685 
1.262 

1.64 

1.45 

1.2 

1.47 

0.67 

1.47 



1.20 
1.24 

1.273 
1.273 



1.45 



QMC 
QMC 

PT a , QMC 
MFT, IPT b 
QMC 

MFT, QMC 
PSCA C 
QMC 
DMFT 

improved Hubbard III 
NRG d 



a Perturbation Theory. 

b Iterated Perturbation Theory. 

c Projective SC Approximation. 

d Numerical Renormalization Group. 
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